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We investigate the effects of scale-free model on cosmology, providing, in this way, a statistical 
background in the framework of general relativity. In order to discuss properties and time evolution 
of some relevant universe dynamical parameters (cosmographic parameters), such as H{t) (Hubble 
parameter), q(t) (deceleration parameter), j{t) (jerk parameter) and a(t) (snap parameter), which 
are well re-defined in the framework of scale-free model, we analyze a comparison between WMAP 
data. Hence the basic purpose of the work is to consider this statistical interpretation of mass 
distribution of universe, in order to have a mass density p dynamics, not inferred from Friedmann 
equations, via scale factor a(t). This model, indeed, has been used also to explain a possible origin 
and a viable explanation of cosmological constant, which assumes a statistical interpretation without 
the presence of extended theories of gravity; hence the problem of dark energy could be revisited in 
the context of a classical probability distribution of mass, which is, in particular, for the scale-free 
model, P{k) ~ k~'^ , with 2 < 7 < 3. The ACDM model becomes, with these considerations, a 
consequence of the particular statistics together with the use of general relativity. 

PACS numbers: 98.80.-k, 98.80.Jk, 98.80.Es 



I. INTRODUCTION 



We found, in literature, many evidences of a lot of un- 
solved problems, which are still open questions in the pic- 
ture of Einstein field equations. One of the most impor- 
tant issue of general relativity (GR) is a complete com- 
prehension and interpretation of the cosmological con- 
stant in accordance with data. In the framework of Fried- 
mann equations we write equations of motion of p, and 
of pressure p, in terms of scale parameter a[t); on the 
other hand we do not consider any kind of mass dynam- 
ics, induced as an emergent effect of a statical mass dis- 
tribution too; this could mean that we lose information 
about matter dynamics and distribution, in sense that 
we consider it as the static part of a dynamic equation, 
where the scale factor a{t) is the true dynamical element. 
In order to consider a more direct matter dynamics we 
focus on p of universe and we supposed the validity of 
the scale-free model into the framework of GR; i.e. to 
say that could exist a statistical model which describes 
the motion of matter and defines implicitly, in this way, 
the true dynamics of the scale factor. This interpretation 
will allow us to describe a dynamics of universe, starting 
only from the basic assumptions that every part of uni- 
verse (in the framework of homogeneity and of isotropy) 
is characterized by a distribution of gravitons, induced 
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by the simple presence of mass. Every graviton is linked 
to each other by a network or by a family of networks 
and its existence is provided by the implicit presence of 
nodes, characterized by a space-time vector , in which 
is possible to have a graviton. This model appears to be 
useful for matter distribution of universe, making us able, 
also, to consider the universe like a disordered system. 

In particular, it is possible to describe the existence 
of heterogeneous spectrum of local interaction patterns 
topologically linked to the presence of nodes. Following 
the features done by other authors Q, we try to inves- 
tigate the consequences of degree heterogeneity on dy- 
namic properties of networks isolating a region of uni- 
verse of correlated networks with a completely random 
organization, defined by the distribution P{k) ~ k^^ ^ 
where 7 has to be considered as an integer number. This 
ansatz shows apparently that dynamics would be limited 
to a particular region but it is generalized, simply, assum- 
ing, as in the beginning, a homogeneous and isotropic 
universe, which allows us to extend a particular result 
to the whole system (the universe evolution). Moreover 
the possibility to consider an elevate number of networks, 
N , is the basic idea to consider firstly, because, although 
< 00, it is true, at the same time, that A^ is a very 
large number. 

In order to describe dynamics, we need all the quan- 
tities [III describing it. These are the cosmographic pa- 
rameters H{t), q{t), j{t) and s{t) P, i, [O, 111, M ill ; 
respectively the Hubble parameter, the deceleration pa- 
rameter, the jerk parameter and finally the snap param- 
eter; these are defined in principle for all t, or redshifts 



2 



z. Assuming the power law degree distribution of the 
form P(fc) ^ ^ which is the basic idea of the model, 
and supposing no interactions among networks, that is 
equivalent to say that the type of interaction described 
by us, i.e. mean-field type, we are also able to infer how 
the scale-free model also deals with transitions and, so, 
it could be also fruitful to explain, for a fixed number of 
gravitons into account, for what values of the involved 
quantities the transition era could start, that could be a 
real topic to test the statistical model proposed [1]. 

It is also possible to characterize the model consider- 
ing the value of parameter 7, which is the most influent 
parameter of the model. For example some evidences de- 
scribe a probability distribution with 2 < 7 < 3, which 
for real systems appears to be a correct choice [6|, 17|. 
Of course, indeed, depending to the value of 7, we find 
many different behaviors; in the present paper we do not 
matter initially about the exact value of 7 in order to 
describe a particular behavior^. 

For our purposes, following the description of 0], 
the first quantity introduced to describe the matter dis- 
tribution is Pfc(t), which represents the average density 
of particles in nodes of degree k and for all the variables, 
it holds 



Pk{t) = -pk + kpzi + A[pk\ 



(1) 



which is a first example of Langevin equation for our 
model; it suggests that the first two terms on the right 
side of equation are due to diffusion in uncorrelated ran- 
dom graphs and the third could be considered as a reac- 
tion kernel that depends only on pk.. From this equation 
we derive the dynamics of a particular choice of k. In 
order to show the evolution of density the density 
which is present in Friedmann equations, we must require 
that 



p(i) = 5]p(fcK(t). 



(2) 



In the continous limit, the above sum becomes an inte- 
grals over the k space, which has to be troncated by using 
k*{t) ~ zi/p{t). Easily, following [1,0] 



P{k)pldk (X 



p^-^ +ckl-^p if kc > k* 
if kc < k* 



k^-'p^ 



(3) 



We consider that near a critical point fic *^=^ k^ den- 
sity is small, thus we expect that the mean-field regime 
holds when k* > kc and N < 00. The corresponding 
Langevin equation is 



p{t) ~ap- bp^ + ^y{ap + 2bp^)/Nr]{t) 



(4) 



where we considered o, 5, c three integration constants 
and, in particular 



kc (X N^/'^. 



(5) 



If we assume universe as a whole system we must require 
the possible presence of noise term; we describe an uni- 
verse without noise if 
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with the uj parameter which has to be in relation with 7 
because, if eq. ^ holds, we must require 



7 < 3w + 3, 



(7) 



or equivalently at least u > —(1/3); because d, 0] w 
appears to be the maximum degree of the network size 
N, it must satisfy a; > 1 and so our interval is always in 
agreement with this request. The basic idea becomes, so, 
to start with the scale-free probability, concerning a dis- 
crete interactions among gravitons and so, after defining 
the density time evolution (eq. (H])), considering the cos- 
mographic evolution of parameters of the universe. All 
the considerations we will do, have to respect eqs. © 
and ([7]). It will be possible to link the universe descrip- 
tion also to critical parameter, as said before, in order 
to describe the transition era, and it will be clear in the 
next sections. 

The article is divided as follows: In the next section we 
describe, briefiy, all the results involved by the use of GR 
in the framework of scale- free model for generic redshift z; 
in particular, we derive the functional form of the dynam- 
ical parameter a(t), first, which will be the basis for the 
construction of the cosmographic parameters, later; their 
own expressions will be found very easily for generic red- 
shifts without approximation in this section. In the third 
section we compare the results, obtained in the previous 
one, in the observable limit z <^ 1 with comparison with 
WMAP data, trying to link also the expressions found 
with the definition of which represents the density of 
cosmological constant A, i.e. dark energy, and we try to 
infer the role of A, just in the statistical way considered. 
The last section is devoted to some qualitative consider- 
ations, conclusion and perspectives for future works on 
this direction. 



II. THE MODEL DYNAMICS 

In this section we infer from scale-free model a specific 
dynamics for universe, without a particular limit on z, or 
of time t, of the dynamical quantities involved. 

Such a dynamics, defined from the statistical back- 
ground is easily valuable, thanks to the first Friedmann 
equation, which is an expression that makes a correspon- 
dence between H{t) and p as follows 



H'^ = ap + P, 



(8) 



See for details 0,[ii|. 



where we have considered, in principle, to take fc = 
in the above expression; and using the positions a = 
and (3 = we can imagine to evaluate p and 

its derivative p by putting them into eq. ([31) , having, 
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so, a differential equation for H{t), in the case r]{t) = 
S{t), which translates the matter statistical distribution 
into an astrophysical differential equation for H(t), which 
reads 



H 



Pf3 bf3^ 



1 



2a H 



2 a J 2a 



f{H) 
N 



(9) 



It could give us the correspondent evolution of universe 
dynamics in terms of H{t) and a{t)] in fact, the solution 
for H{t) is writable, in terms of a{t), because of defini- 
tion H{t) — Jjloga(t). All the discussion is completely 
global, in sense that, at this stage, we are not consider- 
ing any kind of approximation on t, or on z, as already 
mentioned before; that is the same to say that a com- 
plete comprehension of dynamics is, in this paragraph, 
involved. 

For the sake of simplicity we rewrite eq. ((91) at i 7^ 0, 
making the following positions 



a = 



2 2a'' 

p , qP 



b_ 

2a' 



(10) 



and so we get the simple solution for H{t) 



H{t-t*) ^ A^h + Vdt&nh^Vd{t-t*)^, (11) 

where A = the previous expression we consider 

the position d ^ ~ Aac too and we notice that the 
integration constant could be considered as an arbitrary 
time t* . Depending from d sign we could have different 
evolutions of H{t); this sign, indeed, cannot be less than 
zero: A similar choice, i.e. d < 0, shows us a definite time 
interval in which H{t) (and so a{t)) loses its meaning. 
Hence, for describing a variable a{t), for all the time, we 
must have d > 0. On the other hand, this is evident, 
in order to have a solution, which has no periodicity, 
into account. This could, also, be tested, by considering 
directly its expression, d = — Aac and it gets 



P 



-pq(3 



inG 



(12) 



which is always positive for all the possible values as- 
sumed by p and g, if we require that pq > 0. Our choice 
is the simplest one: Because of p, q are " phenomenologi- 
cal" quantities, defining the type of reaction and process, 
we could assume, without losing information, p > and 
q > 0. Indeed, it is possible in principle, that p < 



and q < 0, but we will deny this, comparing also with 
WMAP data. Our efforts, on the other side, describe 
a model, whose validity is strongly dependent from p, q 
values, which has to be chosen by data. 
From integrating eq. ([9]), it is possible to get, also, a 
precise value for H{0). The theoretical expression, that 
we find, is 



1 



j4H^N = p 



ij2 



(3 



2q 



(13) 



which is the solution of equation H{0) — y ^^^^j^^i where 
f{H{t = 0+)) = / (i/(0)), together with H {t = Q-) = 
condition, which is obvious, if we consider the positive 
direction of time. Eq. (|13p gives two positive solutions, 
i.e. 



i7(0)i 



4g/3 

4N-2q 4N-2q 



pa 



-16Npa/3+32Nq(3^ 



4Af-2g 



V2 



H{0)2 = 



pa _ 2ql3 , ^p^a^-16Npal3+32Nqt3^ 
4N-2q 2N-q 2{2N-q) 



V2 



(14) 



(15) 



Both solutions are, in the thermodynamic limit, equiva- 
lent and they would read Hq = 0; in our case, for N < 00, 
only eq. (|15p is relevant, if happens that pa > 2q(3. This 
condition will represent the correct interval for p and q 
choice. 

Moreover, eq. (fT5)) is also useful because tells us the 
value of t* in the expression pT|) ; it reads 



t* = tanh ^ I -=-^ 

yd 



(16) 



After evaluating the dynamics we naturally achieve 
to find the expression for scale factor a(t). We get 



G(t-t*) 

a(t) = a(0)e V2cd ^ 
where the function G {t ~ t*) is expressed by 
G(t-r) - G = 



(17) 



b — \/fitanh ^ 



b + \/dtanh 



lb + y/dta.nh^Vd{t - t*)| 
Vb~^ 

^6-1- \/dtanh|\/d(i - t*)| 
Vb+Vd 



I 



(18) 
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FIG. 1: In this graphic is plotted H{t) with fc = 0. The set 
of parameters is a = 1, 6 = 3 and c — 1, with normalization 
H{0) = 1. 




FIG. 2: In this graphic is plotted a{t) with k = 0. The set 
of parameters is a = 1, & = 3 and c = 1, with normalization 
a(0) = 1. 



with the integration constant G of the form 



G = Vdtanh ^ 



b + \/fitanh 




(19) 



Following all the prescriptions we wrote, and consider- 
ing the dynamical scale factor, in order to describe, for 
all the redshift z, the dynamics of universe, immediately, 
the expression for the deceleration parameter, q{t) could 
be found; it reads 



q{t) = -l- 



c/2d\ b + Vdta.nh Vd(t-t 



[bcosh Vd{t-t*) +\/rfsinh \/d{t-t*) 

(20) 

It appears, of course, negative, which suggests that the 
model is compatible with accelerating universe model. 
The other expressions for j(t) and s(t), which discrim- 
inates the type of universe dynamics, are described in 



FIG. 3: In this graphic is plotted q{t) with k = 0. The set of 
parameters is a = 1, fo = 3 and c — 1. 



appendix for their complexity. 



III. THE OBSERVABLE LIMIT 

Recently the possibility to get, from data, all the ele- 
ments for universe dynamics, has been reported by sev- 
eral works (see for example P, 0, S Hi 20, 21]). A form 
of so-called cosmography, which corresponds to the pos- 
sibility to infer, by some useful parameters, the universe 
evolution, is involved. We know, in fact, that we are 
living in the so-called period of Precision Cosmology, in 
which observations are extremely precise and allows us to 
test better than past, a given evolutionary models. The 
expressions we need are the so-called deceleration q{t), 
jerk j{t), snap s{t) P,[23,[2H. The cosmographic param- 
eters are defined by derivatives of the scale factor for our 
epoch; they, instead of their simplicity, allow to fit the 
data also in terms of redshift, that is an implicit measure 
of time. It is remarkable to remember that the model, 
we are going to consider, here and in the following, is the 
standard model P, d, Q which is not useful when z is not 
little, i.e. z <C 1, that is, indeed, the case of validity of 
the Taylor expansion of scale parameter. 
In fact, for our time, the cosmographic parameters are 
defined by the expression 

^ = l + Ho{t-to)- - t^f + 



and so they reads 



2 • 



H{t) 



lit) 



a 
a 



Jit) 



A3 



(21) 

(22) 
(23) 
(24) 
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sit) 



(25) 



We denote with a subscript these parameters, in sense 
that they are evaluated at t = in; on the meaning of these 
parameters see for example |13| . In this paper we will 
take into account WMAP observations and so we will 
assume that the value of the Hubble constant is Hq ~ 
70 ± 2 km/sec/Mpc at our time t = to ^.\ld^. 

Following these advice, it is clear that, the cosmo- 
CTaphic parameters may also be tested by WMAP data 
[ig] and, above all, for the observed limit z <C 1. It is 
also possible to express them, in terms of the dark en- 
ergy density f^A, which could stress also the role of a 
cosmological constant involved into Einstein field equa- 
tions, and so, give a description of ACDM model in terms 
of statistical background. 

In order to understand such description, we can fol- 
low the prescription of Jl]. Hence, we will use the 
Chevallier - Polarski - Linder (CPL) parameterization for 
the equation of state setting [ij], which reads w — 
wq + Wa — with a{z) = Easily, noting the 

correspondence t < — > z, i.e. dt ^ ~ • Moreover, 

for the scenario, we are considering, in the ACDM case, 
i.e. 



r2,„ + ilk + — 1, 



(26) 



where ftm is normalized^ matter density, flk spatial cur- 
vature density and Ha cosmological constant density, the 
list of above parameters has to be set (wojWa) = (—1,0), 
and, then, after some straightforward calculations, it is 
useful to note qo = ^ — f^^A, jo and sq = 1 — 

Now we have all the ingredients to understand a com- 
plete comparison with scale- free quantities and cosmolog- 
ical parameters. In the next section, so, will be possible 
to write down an expression for A, by comparing with 
the above expression for Q\ in terms of qo . 



So = 1 



2bHo 



6- 



2bHo 



( d 



\2bHQ 



Acd 



(29) 

and inverting, we get expression for the scale-free param- 
eters in function of deceleration, jerk and Hubble param- 
eters as follows 

_ l-so-(l + go)[7 + go + 3(l + go)^] 



, ^ l-.so-(l + go) [7 + go + 3(l + go)^] 



Ml + qoY 



d = -2H^^—^^^ ^ '^"^ [7 + go + 3(1 + qo?] 



4(1 + 9o) 



(32) 



All the previous results deal with a dynamics, induced by 
the use of the scale-free model, and treat the link among 
scale-free parameters and the observed parameters. 
It is necessary to note that all the expressions obtained 
are in the case k = 0, which is the observed one [§]. 
Inducing a curvature term A; ^ is possible to achieve 
an expression for a(t), which is not analytically defined. 
It appears as a strong complication of the model; indeed, 
it seems so also because observation provide a curvature 
density fife ~ 0, but we could easily consider a graphic 
of it, in order to fix ideas. It is important for us to 



a(t) 

4 r 



A. The comparison with WMAP data 

After comparing with the experimental value [1^ [2^ 
defined by Hq, qo,so,la, in a easy way, we link measurable 
parameters, in terms of scale-free parameters 

'» = -'- ik- '''' 

^ Normalized by pc = ■ 



0.0 0.2 0.4 0.6 0.8 1.0 

FIG. 4: In this graphic is plotted a{t) with 7^ 0. The set 
of parameters is a = 1, 6 = 3 and c = 1, with normalization 
a(0) = 1. 

underline that all the graphics are evaluated by a specific 
choice of parameters a, 5, c reported in the graphic labels, 
not chosen by particular reasons but only for graphic 
simplicity. 

Now we can write down all the parameters in terms of 
Hq, as follows 

c - 0.35^0-1 , 

h ~ -2Ho , (33) 
d - 2H?, . 
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and eqs. (?? and [33]) allow us to write down cosmo- 
graphic parameters in terms of Q\ and, in particular, we 
find the definition of A in our model as follows 



cosh^ (^/2i^o^o) 

Remembering the definition of scale-free parameters, 
postulated at the beginning, in particular of critical tran- 
sition parameter, we find for it 

fee 1.429 i?^^. (35) 

In other words the effect of a dark energy are modeled by 
the presence of a cut off fee and for our epoch, in which 
k > kc, all the previous results suggest how to consider 
the correct meaning of A: It could be considered as an 
effect of modeling universe with a P(k) ~ k~^ , or, in 
other words, as a statistical result. 



Moreover, all the previous results, concerning finite 
networks, could be extended in the case of infinite net- 
works by using diverging cut-off kc. It is possible to show, 
that depending from 7 — 1 and q values, different results 
are involved, and this will concern to future efforts. 

Other fields of future interest are to be foimd in the 
possibility to join the statistical behavior, described in 
this work, not necessary in the case of GR 15]. It will 
be clear in future work, that modifications occur if a dif- 
ferent field theory is involved; idea which is possible, be- 
cause the model appears to be independent from the field 
background. Concluding, of course, the statistics in this 
sense is a first kind of example in which the possibility to 
consider a different approach to the problem of A. This 
is possible, also passing through a quantum picture as 
showed, for example, in J^j. The deep difference of this 
kind of work is the classical way in which all the calcula- 
tions have been evaluated. 



IV. CONCLUSION AND PERSPECTIVES 
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We investigated in this paper that the scale- free model, 
achieving a statistical interpretation of the " background" 
of space time, by the presence of a large number of gravi- 
tons, allowed us, thanks to networks interactions, to de- 
fine directly p dynamics, without using, the scale factor 
a(t) dynamics, in the picture of mean field approxima- 
tion. 

Moreover we showed, that, if that dynamics is involved 
into the framework of GR, we are able to define all the 
cosmographic parameters, in terms of time evolution and, 
in principle, also of redshift z. Best results were q{t) < 
which represents the acceleration parameter of universe 
and tells us the universe expansion and also the cosmolog- 
ical constant A, which became dependent from statistical 
and cosmological parameters. This could be thought as 
an explanation of the peculiar reason for what a cosmo- 
logical constant has to be put into Einstein field equa- 
tions, in order to respect the correct statistical matter 
distribution and matter evolution. 

We have, also, discussed the possibility to evaluate the 
parameters of scale-free theory at our epoch in terms of 
measured cosmographic parameters and vice-versa. 

Starting from these results, we can discuss, indeed, also 
the possibility of defining a dark matter dynamical equa- 
tions thanks to the expression of density motion, defined 
in the introduction and, of course, we can also focus a 
possible reformulation of A in the framework of a micro- 
scopic theory, involved by the scale-free model directly. 

The choice of P{k) ~ k~'^ appears to be comparable 
with the density of a high-energy astrophysical object 
and could suggest that scale-free model could be a pos- 
sible adaptable statistical model for the cases of Fermi 
processes [§]. 

The scale-free model appears to be a framework model 
for gravitation so, in principle, extensions of GR are 
adaptable with scale-free model. 
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APPENDIX A: COSMOGRAPHIC PARAMETERS 

The jerk parameter j{t) is defined as follows 



fit) 



Pi + 132 cosh4Vdt + 4 cosh 2\/dt (/?3 -I- 54) 



+(3^ sinh(4Vdt) + 55 + 2Vrfsinh(2Vdt) {fir + gs) : 



cosh"*(\/dt) 
8V2 



b + Vdtanh ( ^/dt 



where f{t) 

g„ = ^liSll T]^e snap parameter s{t) is 
s{t) 



and 



= 511 - .912 + .gi3 + 514 + cosh [&^/dtj (.915 + gie + 517) 



sinh \2^/dtj {(3q + gis - 519 - 320 + 521) 
+ sinh U^fdt] (/3io + .922 — 523 — 524 j 



-I- sinh {&^/dtj (^25 + 525) + 

+2 cosh {A-/dtj {1321 + 528 + 529 - 530 - 53i) + 

- cosh (2\fdtj (/332 - 533 + 534 + 535 + 536). 

The complete list of constants involved into calculations 
is 

[3i = V2(36^ - ihd + Acd?),l32 = \/2h{h^ + 3d), 
/33 = V2{b^ - 2c(f),f3i = 3y/cbd, 13^ = V2d{ib^ + d), 



/36 = I2by/^d, Pr = -V2{-Sb^ + d + Abed), /?« = 6y/^d, 

/dg = -8\/2^d^/2(-962 + 3d+ lObcd), (Al) 
/3io = AV2Zd^^^(9b^ + 3d + 8bcd), = IQb'^, (3i2 = -Ub'^d, 
/3i3 = 2d\Pu = 56bcd^,Pi5 = b\Pie = Qb'd^fin = d^, 
/3i8 = 2053\/d,/3i9 = -126^=^/2/320 = -6462cd3/2^ 
/321 = 56crf^/2 /322 = m^Vd,P23. = -i2b^cA^'^, 



p24 = -~32cd^^^ (325 = ^b^Vd, p26 = ^bd^^^ 
P27 = 2V2cd{3b^ + &(9 + Abc)d + Acd^),[32& = 36^ 

/332 = 16V2cd(-363 + 2^2 cd + Scd^), = -156^ 
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We investigate the effects of scale-free model on cosmology, providing, in this way, a statistical 
background in the framework of general relativity. In order to discuss properties and time evolution 
of some relevant universe dynamical parameters (cosmographic parameters), such as H{t) (Hubble 
parameter), q{t) (deceleration parameter), j(t) (jerk parameter) and s{t) (snap parameter), which 
are well re-defined in the framework of scale-free model, we analyze a comparison between WMAP 
data. Hence the basic purpose of the work is to consider this statistical interpretation of mass 
distribution of universe, in order to have a mass density p dynamics, not inferred from Priedmann 
equations, via scale factor a{t). This model, indeed, has been used also to explain a possible origin 
and a viable explanation of cosmological constant, which assumes a statistical interpretation without 
the presence of extended theories of gravity; hence the problem of dark energy could be revisited in 
the context of a classical probability distribution of mass, which is, in particular, for the scale-free 
model, P{k) ~ fc"'*', with 2 < 7 < 3. The ACDM model becomes, with these considerations, a 
consequence of the particular statistics together with the use of general relativity. 

PACS numbers: 98.80.-k, 98.80.Jk, 98.80.Es 



I. INTRODUCTION 



We found, in literature, many evidences of a lot of un- 
solved problems, which are still open questions in the pic- 
ture of Einstein field equations. One of the most impor- 
tant issue of general relativity (GR) is a complete com- 
prehension and interpretation of the cosmological con- 
stant in accordance with data. In the framework of Pried- 
mann equations wc write equations of motion of p, and 
of pressure p, in terms of scale parameter a{t); on the 
other hand we do not consider any kind of mass dynam- 
ics, induced as an emergent effect of a statical mass dis- 
tribution too; this could mean that we lose information 
about matter dynamics and distribution, in sense that 
we consider it as the static part of a dynamic equation, 
where the scale factor a{t) is the true dynamical clement. 
In order to consider a more direct matter dynamics we 
focus on p of universe and we supposed the validity of 
the scale-free model [7] into the framework of GR; i.e. to 
say that could exist a statistical model which describes 
the motion of matter and defines implicitly, in this way, 
the true dynamics of the scale factor. This interpretation 
will allow us to describe a dynamics of universe, starting 
only from the basic assumptions that every part of uni- 
verse (in the framework of homogeneity and of isotropy) 
is characterized by a distribution of gravitons, induced 
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by the simple presence of mass. Every graviton is linked 
to each other by a network or by a family of networks 
and its existence is provided by the implicit presence of 
nodes, characterized by a space-time vector a;^, in which 
is possible to have a graviton. This model appears to be 
useful for matter distribution of universe, making us able, 
also, to consider the universe like a disordered system. 

In particidar, it is possible to describe the existence 
of heterogeneous spectrum of local interaction patterns 
topologically linked to the presence of nodes. Following 
the features done by other authors [6], we try to inves- 
tigate the consequences of degree heterogeneity on dy- 
namic properties of networks isolating a region of uni- 
verse of correlated networks with a completely random 
organization, defined by the distribution P{k) ~ k~'' , 
where 7 has to be considered as an integer number. This 
ansatz shows apparently that dynamics would be limited 
to a particular region but it is generalized, simply, assum- 
ing, as in the beginning, a homogeneous and isotropic 
universe, which allows us to extend a particular result 
to the whole system (the universe evolution). Moreover 
the possibility to consider an elevate number of networks, 
N, is the basic idea to consider firstly, because, although 
N < 00, it is true, at the same time, that N is a very 
large number. 

In order to describe dynamics, we need all the quan- 
tities [16] describing it. These are the cosmographic pa- 
rameters H{t), q{t), jit) and s{t) [1, 9, 17, 18, 20, 21]; 
respectively the Hubble parameter, the deceleration pa- 
rameter, the jerk parameter and finally the snap param- 
eter; these are defined in principle for all t, or redshifts 



z. Assuming the power law degree distribution of tlie 
form P(fc) ~ A;"''', which is the basic idea of the model, 
and supposing no interactions among networks, that is 
equivalent to say that the type of interaction described 
by us, i.e. mean-field type, we are also able to infer how 
the scale-free model also deals with transitions and, so, 
it could be also fruitful to explain, for a fixed number of 
gravitons into account, for what values of the involved 
quantities the transition era could start, that could be a 
real topic to test the statistical model proposed [5]. 

It is also possible to characterize the model consider- 
ing the value of parameter 7, which is the most influent 
parameter of the model. For example some evidences de- 
scribe a probability distribution with 2 < 7 < 3, which 
for real systems appears to be a correct choice [6, 7]. 
Of course, indeed, depending to the value of 7, we find 
many different behaviors; in the present paper we do not 
matter initially about the exact value of 7 in order to 
describe a particular behavior^. 

For our purposes, following the description of [6, 7], 
the first quantity introduced to describe the matter dis- 
tribution is Pk{t), which represents the average density 
of particles in nodes of degree k and for all the variables, 
it holds 



Pk{t) = -Pk + kpzi + A[pk], 



(1) 



which is a first example of Langevin equation for our 
model; it suggests that the first two terms on the right 
side of equation are due to diffusion in uncorrelated ran- 
dom graphs and the third could be considered as a reac- 
tion kernel that depends only on pk- From this equation 
we derive the dynamics of a particular choice of k. In 
order to show the evolution of density p{t), the density 
which is present in Priedmann equations, we must require 
that 



p{t) = ^P{k)pk{t). 



(2) 



In the continous limit, the above sum becomes an inte- 
grals over the k space, which has to be troncated by using 
k*{t) ~ zi/p{t). Easily, following [6, 7] 



I P{k)pldk (X I 



pT-i + ckl-'p 



if kc > k* 
if kr < k* 



(3) 



We consider that near a critical point pc kc, den- 
sity is small, thus we expect that the mean-field regime 
holds when k* > kc and N < 00. The corresponding 
Langevin equation is 



p{t) c^ap- bp" + ^{ap + 2bp-^)/NTi{t) (4) 

where we considered a, b, c three integration constants 
and, in particular 



kcOcN^^"^. 



(5) 



If we assume universe as a whole system we must require 
the possible presence of noise term; we describe an uni- 
verse without noise if 



(6) 



with the u! parameter which has to be in relation with 7 
because, if eq. (6) holds, we must require 



7 < 3w + 3, 



(7) 



or equivalently at least u) > —(1/3); because [6, 7] to 
appears to be the maximum degree of the network size 
N, it must satisfy uj > 1 and so our interval is always in 
agreement with this request. The basic idea becomes, so, 
to start with the scale-free probability, concerning a dis- 
crete interactions among gravitons and so, after defining 
the density time evolution (eq. (4)), considering the cos- 
mographic evolution of parameters of the universe. All 
the considerations we will do, have to respect eqs. (6) 
and (7). It will be possible to link the universe descrip- 
tion also to critical parameter, as said before, in order 
to describe the transition era, and it will be clear in the 
next sections. 

The article is divided as follows: In the next section we 
describe, briefly, all the results involved by the use of GR 
in the framework of scale- free model for generic redshift z; 
in particular, we derive the functional form of the dynam- 
ical parameter a{t), first, which will be the basis for the 
construction of the cosmographic parameters, later; their 
own expressions will be found very easily for generic red- 
shifts without approximation in this section. In the third 
section we compare the results, obtained in the previous 
one, in the observable limit z <C 1 with comparison with 
WMAP data, trying to link also the expressions found 
with the definition of Oa, whic;li represents the density of 
cosmological constant A, i.e. dark energy, and we try to 
infer the role of A, just in the statistical way considered. 
The last section is devoted to some qualitative consider- 
ations, conclusion and perspectives for future works on 
this direction. 



II. THE MODEL DYNAMICS 

In this section we infer from scale-free model a specific 
dynamics for universe, without a particular limit on z, or 
of time t, of the dynamical quantities involved. 

Such a dynamics, defined from the statistical back- 
ground is easily valuable, thanks to the first Friedmann 
equation, which is an expression that makes a correspon- 
dence between H{t) and p as follows 



ap + (3, 



(8) 



1 See for details [7, 19]. 



where we have considered, in principle, to take A; = 
in the above expression; and using the positions a = 
and P = we can imagine to evaluate p and 

its derivative p by putting them into eq. (4), having, 



so, a differential equation for H{t), in the case ri{t) = 
5{t), which translates the matter statistical distribution 
into an astrophysical differential equation for H{t), which 
reads 

'pf3 b(3'\ 1 



+ 



q(3 



H 



2a 



f{H) 
N 



5{t). 



(9) 



It could give us the correspondent evolution of universe 
dynamics in terms of H{t) and a(t); in fact, the solution 
for H{t) is writable, in terms of a(t), because of defini- 
tion H{t) — -^\oga{t). All the discussion is completely 
global, in sense that, at this stage, we are not consider- 
ing any kind of approximation on or on 0, as already 
mentioned before; that is the same to say that a com- 
plete comprehension of dynamics is, in this paragraph, 
involved. 

For the sake of simplicity we rewrite eq. (9) at t 7^ 0, 
making the following positions 



a = 



p(3 



'20.' 



2 a 



(10) 



b 

and so we get the simple solution for H{t) 



H{t-t*)=A^b + Vdta,nh^^Vd{t-t*)Y (11) 

where A = In the previous expression we consider 

the position d = b"^ — 4ac too and we notice that the 

integration constant could be considered as an arbitrary 
time t* . Depending from d sign we could have different 
evolutions of H{t); this sign, indeed, cannot be less than 
zero: A similar choice, i.e. d < 0, shows us a definite time 
interval in which H{t) (and so a{t)) loses its meaning. 
Hence, for describing a variable a{t), for all the time, we 
must have d > 0. On the other hand, this is evident, 
in order to have a solution, which has no periodicity, 
into account. This could, also, be tested, by considering 
directly its expression, d = b^ — Aac and it gets 



(12) 



which is always positive for all the possible values as- 
sumed by p and q, if we require that pq > 0. Our choice 

is the simplest one: Because of p, q are " phenomenologi- 
cal" quantities, defining the type of reaction and process, 
we could assume, without losing information, j; > and 
g' > 0. Indeed, it is possible in principle, that p < Q 



and q < 0, but we will deny this, comparing also with 
WMAP data. Our efforts, on the other side, describe 
a model, whose validity is strongly dependent from p, q 
values, which has to be chosen by data. 
From integrating eq. (9), it is possible to get, also, a 
precise value for H{0). The theoretical expression, that 
we find, is 



a 



f(H(0)) 



N ' 



(13) 



where 



which is the solution of equation H(0) - 

f{H{t = 0+)) = / (ff (0)), together with F (t = 0") = 
condition, which is obvious, if we consider the positive 
direction of time. Eq. (13) gives two positive solutions, 
i.e. 



if(0)i = 



pa 



4N-2q 



^■p^cx^ - 16Afpa/3+32iVg/32 
4N-2q 4JV-2q 



V2 



pa _ 2ql3 , ^/p^a^ -16Npat3+32Nqf3^ 
4N-2q 2N-q 2{2N-q) 



V2 



(14) 



.(15) 



Both solutions are, in the thermodynamic limit, equiva- 
lent and they would read i?o = 0; in our case, for N < 00, 
only eq. (15) is relevant, if happens that pa > 2q(3. This 
condition will represent the correct interval for p and q 
choice. 

Moreover, eq. (15) is also useful because tells us the 
value of t* in the expression (11); it reads 



r=tanh- ('-^ 



(16) 



After evaluating the H{t) dynamics we naturally achieve 
to find the expression for scale factor a{t). We get 



G({-t*) 

a{t) = a(0)e V2c<i , 
where the function G [t — t*) is expressed by 
G{t-t*)-G = 



(17) 



I 4/6-t-\/rftanh|\/d(f-t*)| 

= -yb- v^tanh"^ < > 

Vb-Vd 



+ 



+ y b + Vdtanh ^ < 



b + vdtanh (i - t*)| 
Vb+Vd 



> (18) 
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FIG. 1: In this graphic is plotted H{t) with A; = 0. The set 
of parameters is o = 1, 6 = 3 and c = 1, with normalization 
H{0) = 1. 




1.0 



FIG. 2: In this graphic is plotted a{i) with fc = 0. The set 
of parameters is o = 1, 6 = 3 and c = 1, with normalization 
a(0) = 1. 



with the integration constant G of the form 



G = Vdtanh ^ 



-\Jh + \/rftanh 



Vb+Vd 



(19) 



Following all the prescriptions we wrote, and consider- 
ing the dynamical scale factor, in order to describe, for 
all the redshift z, the dynamics of universe, immediately, 
the expression for the deceleration parameter, q{t) could 
be found; it reads 



lit) 



s/c/2d^b+Vdtanh Vd{t-t* 



(^6 cosh Vd{t-t*) +V^sinh Vd{t-t*) 

(20) 

It appears, of course, negative, which suggests that the 
model is compatible with accelerating universe model. 
The other expressions for j(t) and s(t), which discrim- 
inates the type of universe dynamics, are described in 



FIG. 3: In this graphic is plotted q{t) with k = 0. The set of 
parameters is o = 1, 6 = 3 and c = 1. 



appendix for their complexity. 



III. THE OBSERVABLE LIMIT 

Recently the possibility to get, from data, all the ele- 
ments for universe dynamics, has been reported by sev- 
eral works (see for example [1-4, 20, 21]). A form of 
so-called cosmography, which corresponds to the possi- 
bility to infer, by some useful parameters, the universe 
evolution, is involved. We know, in fact, that wc are 
living in the so-called period of Precision Cosm,ology, in 
which observations are extremely precise and allows us to 
test better than past, a given evolutionary models. The 
expressions we need are the so-called deceleration q(i). 
jerk j(t), snap s{t) [1, 20, 21]. The cosmographic param- 
eters are defined by derivatives of the scale factor for our 
epoch; they, instead of their simplicity, allow to fit the 
data also in terms of redshift, that is an implicit measure 
of time. It is remarkable to remember that the model, 
we are going to consider, here and in the following, is the 
standard model [1, 8, 9] which is not useful when z is not 
little, i.e. 2; <C 1, that is, indeed, the case of validity of 
the Taylor expansion of scale parameter. 
In fact, for our time, the cosmographic parameters are 
defined by the expression 

^ = 1 + i/o(t - to) - U^Hl{t - tof + 
ao I 



and so they reads 



m = 



2 • 



jit) 



(21) 

(22) 

(23) 
(24) 



s{t) 



(25) 



Wc denote with a subscript these parameters, in sense 
that they are evaluated a.t t = to; on. the meaning of these 
parameters see for example [13]. In this paper we will 
take into account WMAP observations and so we will 
assume that the value of the Hubble constant is Hq ~ 
70 ± 2 km/sec/Mpc at our time t = to [9, 10]. 

Following these advice, it is clear that, the cosmo- 
graphic parameters may also be tested by WMAP data 
[10] and, above all, for the observed limit z ^ 1. It is 
also possible to express them, in terms of the dark en- 
ergy density flA, which could stress also the role of a 
cosmological constant involved into Einstein field equa- 
tions, and so, give a description of ACDM model in terms 
of statistical background. 

In order to understand such description, we can fol- 
low the prescription of [11]. Hence, we will use the 
Chevallier-Polarski-Linder (CPL) parameterization for 
the equation of state setting [14], which reads w = 
wq + Wa — jq^, with a{z) = ji^- Easily, noting the 

correspondence t < — > z, i.e. dt = — ^^^^^y^- Moreover, 
for the scenario, we are considering, in the ACDM case, 
i.e. 



(26) 



where flm is normalized^ matter density, Qk spatial cur- 
vature density and Cl\ cosmological constant density, the 
list of above parameters has to be set {wa,Wa) = (—1,0), 
and, then, after some straightforward calculations, it is 
useful to note qo = ^ — |^^a, jo and sq = 1 — |fim. 

Now we have all the ingredients to understand a com- 
plete comparison with scale- free qiiantities and cosmolog- 
ical parameters. In the next section, so, will be possible 
to write down an expression for A, by comparing with 
the above expression for Oa in terms of qo- 



So 



Acd 

2bHoJ r ~ \2bHo ) ^ [mi^ ) ~ "T 

(29) 

and inverting, we get expression for the scale-free param- 
eters in function of deceleration, jerk and Hubble param- 
eters as follows 



c = 



1 - So - (1 + qo) [7 + qo+ 3(1 + qof 



8Hoil + qof 



b = Ho 



l-SQ- (1 + go) [7 + go + 3(l + go)^] 
4(l + go)2 



l-so-(l + go) [7 + go + 3(l + go)^ 
4(1 + qo) 



(30) 



(31) 



(32) 



All the previous results deal with a dynamics, induced by 
the use of the scale-free model, and treat the link among 
scale-free parameters and the observed parameters. 
It is necessary to note that all the expressions obtained 
are in the case k = 0, which is the observed one [9]. 
Inducing a curvature term k ^ is possible to achieve 
an expression for a{t), which is not analytically defined. 
It appears as a strong complication of the model; indeed, 
it seems so also because observation provide a curvature 
density fife ~ 0, but we could easily consider a graphic 
of it, in order to fix ideas. It is important for us to 




A. The compEirison with WMAP data 

After comparing with the experimental value [23, 24] 
defined by Hq, qa, sq, Iq, in a easy way, we link measurable 
parameters, in terms of scale-free parameters 

Normalized by pc = 



0.0 0.2 0.4 0.6 0.8 1.0 

FIG. 4: In this graphic is plotted a(t) with fc 0. The set 
of parameters is o = 1, 6 = 3 and c = 1, with normalization 
a(0) = 1. 

underline that all the graphics are evaluated by a specific 
choice of parameters o, 6, c reported in the graphic labels, 
not chosen by particular reasons but only for graphic 
simplicity. 

Now we can write down all the parameters in terms of 
i?0) as follows 

c ~ 0.35iIo"^ , 

h ~ -2i?o , (33) 
d ~ 2i?Q . 



and eqs. (?? and 33) allow us to write down cosmo- 
graphic parameters in terms of and, in particular, we 
find the definition of A in our model as follows 



cosh^ (\/2ffoio) ' 

Remembering the definition of scale- free parameters, 
postulated at the beginning, in particular of critical tran- 
sition parameter, we find for it 

fee ~ 1.429 if^. (35) 

In other words the effect of a dark energy arc modeled by 
the presence of a cut off kc and for our epoch, in which 
k > kc, all the previous results suggest how to consider 
the correct meaning of A: It could be considered as an 
effect of modeling universe with a P{k) ~ k~'^ , or, in 
other words, as a statistical result. 



Moreover, all the previous results, concerning finite 
networks, could be extended in the case of infinite net- 
works by using diverging cut-off kc- It is possible to show, 
that depending from 7 — 1 and q values, different results 
are involved, and this will concern to future efforts. 

Other fields of future interest are to be found in the 
possibility to join the statistical behavior, described in 
this work, not necessary in the case of GR [15]. It will 
be clear in future work, that modifications occur if a dif- 
ferent field theory is involved; idea which is possible, be- 
cause the model appears to be independent from the field 
background. Concluding, of course, the statistics in this 
sense is a first kind of example in which the possibility to 
consider a different approach to the problem of A. This 
is possible, also passing through a quantum picture as 
showed, for example, in [22]. The deep difference of this 
kind of work is the classical way in which all the calcula- 
tions have been evaluated. 



IV. CONCLUSION AND PERSPECTIVES 
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We investigated in this paper that the scale- free model, 
achieving a statistical interpretation of the " background" 
of space time, by the presence of a large number of gravi- 
tons, allowed us, thanks to networks interactions, to de- 
fine directly p dynamics, without using, the scale factor 
a{t) dynamics, in the picture of mean field approxima- 
tion. 

Moreover we showed, that, if that dynamics is involved 
into the framework of GR, we are able to define all the 
cosmographic parameters, in terms of time evolution and, 
in principle, also of redshift z. Best results were q{t) < 
which represents the acceleration parameter of universe 
and tells us the universe expansion and also the cosmolog- 
ical constant A, which became dependent from statistical 
and cosmological parameters. This could be thought as 
an explanation of the peculiar reason for what a cosmo- 
logical constant has to be put into Einstein field equa- 
tions, in order to respect the correct statistical matter 
distribution and matter evolution. 

We have, also, discussed the possibility to evaluate the 
parameters of scale-free theory at our epoch in terms of 
measured cosmographic parameters and vice-versa. 

Starting from these results, we can discuss, indeed, also 
the possibility of defining a dark matter dynamical equa- 
tions thanks to the expression of density motion, defined 
in the introduction and, of course, we can also focus a 
possible reformulation of A in the framework of a micro- 
scopic theory, involved by the scale-free model directly. 

The choice of P{k) ~ k~'^ appears to be comparable 
with the density of a high-energy astrophysical object 
and could suggest that scale-free model could be a pos- 
sible adaptable statistical model for the cases of Fermi 
processes [9]. 

The scale-free model appears to be a framework model 
for gravitation so, in principle, extensions of GR are 
adaptable with scale-free model. 
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APPENDIX A: COSMOGRAPHIC PARAMETERS 

The jerk parameter j{t) is defined as follows 

m = A /32 cosh4\/rft + 4 cosh2\/dt (/?3 + 54) + 

+(35 smh{4:Vdt) +96 + 2Vrfsinh(2Vrft) (/Jy + ffs) , 



cosh-''(v^t) 



where f{t) 

Qn = ^2^1*1. The snap parameter s{t) is 



b + Vdtanh (Vdt^ 



and 



-I- 



= 911 - 912 + 913 + 514 + cosh ^eVdt^ {gi5 + gie + gn) 
sinh {2Vdtj {Pg + gis - gw - 520 + 521) + 
+ sinh (AVdt^ (/3io -|- 522 - 523 - 524) 

-I- sinh (eVdtj (525 + 526) + 

+2 cosh (^AVdt^ {1327 + 928 + 529 - 530 - 531) 

- cosh (2\/dtj (/332 - 533 + 534 + 535 + 536). 

The complete list of constants involved into calculations 
is 

/Ji = \/2(36^ - 3bd + 4crf2), /J2 = V2b{b'^ + 3d), 
03 = V2{b^ - 2c(f), I3i = 3y/cbd, 05 = ^2^(36^ + d), 



+ 



Pe = 12by/cd, 07 = -\/2(-36^ + d + Abed), Ps = 6\/crf, 

/Jg = -8\/2crf^/^(-96^ + 3d+ Wbcd), (Al) 
Pio = 4V2cd3/2(-g^2 + + gijcd), = 106'', /3i2 = -1262^, 
/3i3 = 2d2,/3i4 = 56bcd^,Pi5 = b\Pw = 6b^d,f3i7 = d^, 
/3i8 = 20b^Vd,l3i9 = -126d^/^/J2o = -646^ cc^^/^ 
/321 = 56cc^5/^/322 = 1663\/rf,/323 = -3262cc^^/^ 



/324 = -32crf^/^/?25 = Ab^Vd, p26 = Abd^''^, 
p27 = 2V2cd{3b^ + 6(9 + 46c)rf + Acd"^), /J28 = 36*, 

P29 = 662(i,/330 = -rf^/?3l = -326crf^ 
/?32 = 16\/2cd(-36^ + 26^crf + 3c(i^),/333 = -156*, 
/334 = 662rf,/335 = d^/336 = 86crf2. 
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